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The distribution of dark matter in the Galaxy, according to state-of-the-art simulations, shows 
not only a smooth halo component but also a rich substructure where a hierarchy of dark matter 
subhalos of different masses is found. We present a search for potential dark matter subhalos 
in our Galaxy exploiting the high (HE, 100 MeV - 100 GeV) and very-high-energy (VHE, >100 
GeV) y-ray bands. We assume a scenario where the dark matter is composed of weakly interacting 
massive particles of mass over 100 GeV, and is capable of self-annihilation into standard model 
products. Under such a hypothesis, most of the photons created by the annihilation of dark matter 
particles are predicted to lay in the HE y-ray band, where the Fermi-Large Area Telescope is 
the most sensitive instrument to date. However, the distinctive spectral cut-off located at the 
dark matter particle mass is expected in the VHE y-ray band, thus making imaging atmospheric 
Cherenkov telescopes like VERITAS the best suited instruments for follow-up observations and 
the characterization of a potential dark matter signature. We report on the ongoing VERITAS 
program to hunt for these dark matter subhalos, particularly focusing on two promising dark 
matter subhalo candidates selected among the Fermi-LAT Second Source Catalog unassociated 
high-energy y-ray sources. 
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1. Introduction 

There are strong observational evidences that support the existence of dark matter (DM) in 
our universe, like the dynamics of galaxies and galaxy clusters, and the gravitational lensing effect. 
The abundance of this unknown component of the universe has been quantified to 84% of its total 
mass-density [1], Assuming the hypothesis of DM being composed of weakly interacting massive 
particles (WIMPs), annihilation or decay of DM into Standard Model particles opens the possibility 
of indirectly detecting DM through the observation of these products, amongst them y-ray photons 
(see [2]). The expected flux of prompt y-ray emission from DM annihilation features a distinctive 
spectral shape primarily characterized by a cut-off at the WIMP mass. The preferred targets to 
search this DM spectral signature from are those regions showing the highest signal strength, by 
holding a large concentration of DM and/or being located close to the observer, like the Galactic 
Center and Halo, the dwarf spheroidal satellites of the Milky Way, and galaxy clusters. No clear 
DM signal has been detected so far in any of those objects [2]. 

High-resolution simulations of Milky Way-like DM halos indicate the presence of substruc¬ 
ture down to the smallest resolved scales (see, e.g., [3]). Not all of these subhalos may have ac¬ 
cumulated significant amounts of baryonic matter and would therefore be invisible to astronomical 
observations from radio to X-ray energies. Assuming a self-annihilating WIMP in the GeV-TeV 
mass range, these subhalos would be only visible at y-ray energies. Since the emission from WIMP 
annihilation is expected to be steady, such hypothetical sources would be found in deep sky y-ray 
surveys, and could already be part of the collection of Fermi -Large Area Telescope (LAT) sources 
showing no conventional counterpart at any other wavelengths [4]. If the distinct spectral cut-off 
at the WIMP mass is located at energies too high to be measurable by Fermi-L/VT, imaging at¬ 
mospheric Cherenkov telescopes (IACTs) like VERITAS may be crucial to identify such a DM 
signature. Several searches for DM subhalos in the HE y-ray band have been conducted (see, e.g., 
[5, 6, 7]) and some of them have already triggered follow-up observations by IACTs [8, 9]. 

This contribution focuses on the observations of two DM subhalo candidates in the HE and 
VHE y-ray bands. Section 2 describes the source selection, while the VERITAS VHE y-ray and 
Fermi -LAT HE y-ray observations and data analyses are contained in Sections 3 and 4 respectively. 
In Section 5 the results are discussed. A brief summary and outlook can be found in Section 6. 

2. Source selection 

The Second Fermi -LAT Catalog (2FGL) contains 1873 HE y-ray sources detected by the LAT 
instrument after the first 24 months of observations, 576 sources lacking any clear association. 
We adapted the selection criteria from [8] to find good DM subhalo candidates among the latter 
unassociated sources. We requested our candidates: I. to be located at high Galactic latitudes (\b\ > 
10°); II. not to present a variable flux; III. not to have potential counterparts at other wavelengths; 
IV. to be observable from VERITAS latitude with a maximum culmination zenith angle of 40°. 
To satisfy criterium III. we examined the 95% containment region of the candidates, looking for 
cataloged sources in the HEASARC Archive 1 and analyzed Swift-X RT 2 exposures when available. 
In addition, we estimated the required observation time for a 5a detection with VERITAS, based 

1 http://heasarc.gsfc.nasa.gov/ 2 http://swift.gsfc.nasa.gov/ 
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2FGL J1115.0-0701 (right panel). Both significance distributions are well described by the hypothesis of 
background fluctuations, as illustrated by the significance histograms insets. 

on a 2FGL Catalog flux extrapolation to the VHE range, discarding candidates that would require 
more than 50 hours of VERITAS time to be detected. We ended up with two best candidates, 
namely, 2FGL J0545.6+6018 and 2FGL J1115.0—0701. It is worth stressing that, by the time 
these sources were proposed for VERITAS observations, the 2FGL Catalog was the most updated 
publicly released catalog of Fermi -LAT point like sources. Thus, 2FGL Catalog denominations for 
our DM subhalo candidates will be used throughout the text. 

3. VERITAS observations and data analysis 

The Very Energetic Radiation Imaging Telescope Array System [10] is a ground-based y-ray 
telescope array located at the Fred Lawrence Whipple Observatory in southern Arizona (31 40N, 
110 57W, 1.3 km a.s.l.). The array consists of four imaging atmospheric-Cherenkov telescopes, 
each employing a tessellated 12 m Davies-Cotton reflector instrumented with a photomultiplier- 
tube camera with a 3.5° field of view. VERITAS is designed to detect emission from astrophysical 
objects in the energy range from 85 GeV to greater than 30 TeV, with a nominal sensitivity suf¬ 
ficient to detect, at the 5 a level, a steady point-like source with 1% of the Crab Nebula flux in 
approximately 25 hrs. VERITAS has an energy resolution of 15% and an angular resolution (68% 
containment) of 0.1° at 1 TeV. 

The VERITAS observations of 2FGL J0545.6+6018 and 2FGL J1115.0—0701 presented here 
were made under clear, moonless skies between October 2013 and March 2015, all of them after 
the 2012 upgrade of the VERITAS cameras and trigger system. After data-quality selection, 2FGL 
J0545.6+6018 observations amount to 8.5 hrs performed at an average elevation of 61°, while 2FGL 
J1115.0—0701 observations span 13.8 hrs at an average elevation of 51°. The data were acquired in 
the so called wobble (also known as false source ) observation mode. The shower images from each 
telescope are calibrated, cleaned, and parametrized using a momentum analysis. These parameters 
are utilized to discriminate y-ray like events against hadronic events applying multivariate analysis 
techniques, and are also combined to estimate the arrival direction and energy of the original y-ray. 
The event selection cuts are optimized a priori over a Crab Nebula data set. 
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Source 

E, h [GeVl 

N on 

N off 

Nexc 

G 

m UL{ 95% c.l.) 

^UL{ 99 % c.l.) 


150 

327 

5836 

-16.2 ± 18.6 

-0.9 

24.0 

44.1 

2FGL J0545.6+6018 

250 

94 

1704 

-6.7±10.0 

-0.7 

28.4 

52.0 


450 

24 

536 

-8.1±5.1 

-1.5 

3.3 

8.4 


150 

673 

10024 

-4.7 ± 26.8 

-0.2 

99.7 

59.2 

2FGL J1115.±0-0701 

250 

144 

2112 

0.5 ± 12.4 

0.0 

53.6 

32.0 


450 

31 

488 

-2.6 ± 5.8 

-0.4 

21.0 

12.2 


Table 1: Upper limit to the number of excess events for three different energy thresholds £)/,. N on and 
N a ff are the number of gamma-like events in the signal and background regions respectively, N exc being the 
number of excess event in the signal region. The significance a is computed applying Eq. 17 in [12], 






Integral 




Differential 



E>150 GeV 

E>250 GeV 

E>450 GeV 

E 0 = 1 TeV 


c.L 

[cm 2 s *] 

[C.U.] 

[cm 2 s '] 

[C.U.] 

[cm 2 s '] 

[C.U.] 

[TeW'cm^s -1 ] 

2FGL J0545.6+6018 

95% 

1.95xl0~ 12 

3.57xl0~ 12 

0.6% 

0.95x 10~ 12 
1.69x 10~ 12 

0.6% 

0.16x 10~ 12 
0.42x 10~ 12 

0.2% 

1.88xl0- 13 
2.38x 10~ 13 

99% 

1.0% 

1.0% 

0.6% 

2FGLJ1115.0-0701 

95% 

2.15xl0~ 12 

3.62xl0~ 12 

0.6% 

0.88x 10~ 12 
1.48x 10~ 12 

0.5% 

0.31x 10~ 12 
0.53xl0- 12 

0.4% 

2.11x 10~ 13 
3.54x 10~ 13 

99% 

1.0% 

0.9% 

0.8% 


Table 2: Upper limits to the integral fluxes for three different energy thresholds. The upper limits to the 
differential flux are computed for a pivot energy of 1 TeV. 


We find no evidence for y-ray emission above VERITAS energy threshold in any of the ob¬ 
served fields, as illustrated by the significance maps shown in Fig. 1. We note that the distribution 
of significances in either fields are compatible with pure background fluctuations. It is worth men¬ 
tioning that the sources’ coordinates in the 3FGL Catalog do not significantly change with respect 
to the values provided in the 2FGL Catalog (and used for VERITAS observations). The angular 
distance between both sets of coordinates are 3 arcmin and 4 arcsec for 2FGE J0545.6+6018 and 
2FGE J1115.0—0701 respectively, well below VERITAS angular resolution. 

The results from the VERITAS data analyses are summarized in Table 1 , where observations 
significances and upper limits to the number of excess events coming from each source can be 
found, and Table 2, where integral and differential upper limits to the VHE y-ray flux coming 
are listed. For the calculation of the upper limits to the number of excess events we applied the 
Rolke et al. method [11], For the calculation of the upper limits to the flux we assumed power-law 
spectra with spectral indices of -2.0 and -2.1 for 2FGL J0545.6+6018 and 2FGL J1115.0—0701 
respectively 3 . We constrain the integral flux of the sources in the VHE y-ray domain to be below 
1% of the Crab Nebula integral flux at a 95% confidence level. 


4. Fermi -LAT observations and data analysis 

The LAT is a space-based electron-positron pair-conversion instrument on board the Fermi 
y-ray Space Telescope [13]. It is sensitive to yrays from 20 MeV to more than 300 GeV. The LAT 
has a field of view of 2.4 sr and operates primarily in an all-sky survey mode, covering the entire 
sky approximately every three hours. 

We examined 6.7 years of LAT data (2008-07-31 to 2015-04-02, dubbed 7-years analysis 
henceforth) selecting reprocessed Pass 7 events in the energy range between 100 MeV and 300 GeV. 

3 These spectral indices are the best fit values to power law spectra from dedicated LAT analyses presented in Sec. 4. 
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Analysis 

E 0 

N 

a 

P 



[MeV] 

[ x 10~ 13 MeV -1 cm~ 2 s“‘; 

i 



2FGL 

2292 

(1.6 ± 0.3)x 10~ 13 

1.95 + 0.11 

n/a 

2FGL J0545.6+6018 

3FGL 

1913 

(4.0 ± 0.5)x 10~ 13 

1.79 + 0.13 

0.30 ± 0.09 


7-Years 

1913 

(3.7 ± 0.4)x 10~ 13 

1.76 + 0.12 

0.33 ± 0.08 


2FGL 

7491 

(7.8 ± 2.7)x 10~ 15 

1.60 + 0.28 

n/a 

2FGLJ1115.0-0701 

3FGL 

1947 

(1.7 ± 0.2)x 10~ 13 

2.11 ±0.11 

n/a 


7-Years 

1947 

(1.3 ± 0.1)x 10~ 13 

2.13 + 0.10 

n/a 


Table 3: Spectral parameters for 2FGL J0545.6+6018 and 2FGL J1115.0—0701 coming from the 2FGL and 
3FGL Catalogs as well as our dedicated 7-year analysis. For sources showing a power-law best fit to their 
spectra these are parametrized as dN/dE = N(E /Ef while for sources best described by a log-parabola 
their spectra are parametrized as dN/dE = N(E/Ef)^ a ^P ,og ^ E / Eo \ All errors are statistical. 
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Figure 2: Fermi-LAT spectrum and 1 cr confidence contour of the fit, together with the VERITAS differential 
flux upper limits (95% c.l.) for 2FGL J0545.6+6018 (left) and 2FGL J1115.0—0701 (right). The Crab 
Nebula differential spectrum in the VHE y-ray band [14] at 100%, 10%, and 1% levels is shown as reference. 


We selected class 2 events recorded at zenith angles below 100° during good time intervals. We 
analyzed regions of interest of 10° radius centered on the updated 3FGL coordinates of the sources. 
The data were analyzed using the Fenni Science Tools version v9r33p0-fssc-2014 0520. 

Our 7-year analysis confirms the spectral characterization found in the 3FGL Catalog for the 
two sources of interest, and disfavors the spectral characterization found in the 2FGL Catalog. The 
collection of spectral parameters for both 2FGL J0545.6+6018 and 2FGL J1115.0—0701 are shown 
in Table 3. The agreement between our analyses and the 3FGL Catalog is also illustrated in Fig. 2. 

In addition, we independently confirm the variable nature of the y-ray emission from 2FGL 
J1115.0—0701 as first reported in the 3FGL Catalog. More specifically, the source presented a 
high activity state in mid January 2011 lasting for approximately two weeks with integral fluxes of 
~ 4 x 10~ 9 cm~ 2 s -1 above 1 GeV (TS = 37), while its TS value for the last 12 months of analyzed 
data is below 10 for energies larger than 100 MeV. The difference (> 1 a) in the normalization factors 
between our analysis and the 3FGL Catalog nominal value could be attributed to the apparent fade 
off of 2FGL J1115.0—0701, while the spectral indices are compatible within errors. 


5. Results 

The flux variability for 2FGL J1115.0—0701 first noted in the 3FGL Catalog and later verified 
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by our dedicated LAT analysis completely rejects the hypothesis of this source being a DM subhalo. 
Our 7-years analysis showed a high flux state with approximately two weeks duration compatible 
with an integral flux above 1 GeV of ~ 4 x 10 9 cm~ 2 s” 1 (14% of the inverse Compton component 
of the Crab Nebula in the HE y-ray band). VERITAS constraints to 2FGL J1115.0—0701 VHE 
flux would rule out the strict extrapolation of the 2FGL and 3FGL spectra into the VHE domain 
as described in Table 3. The extrapolation of the 7-years spectral characterization is not ruled 
out by the VERITAS constraints but it is in strong tension with them. In addition, VERITAS upper 
limits are still compatible within the extrapolated 1 a confidence contour coming from each of those 
analyses. The tension between the VERITAS limits and our 7-years LAT analysis might favor some 
spectral softening at VHE which, if attributed to attenuation due to extragalactic background light 
absorption, together with the variable HE flux of the source and its high Galactic latitude, might 
support a distant blazar hypothesis as the nature of this source. However, it is also possible that 
the source was actually in a low state during part or all of the VERITAS observations, correlated 
to the lack of activity observed during the last 12 months of analyzed LAT data. In such a scenario 
nothing could be said in regards to the spectral softening of 2FGL J1115.0—0701 in the VHE band. 

While VERITAS limits rule out the direct extrapolation to the VHE domain of the 2FGL 
Catalog spectral description of 2FGL J0545.6+6018, they do not constrain at ah the log parabola 
description from the 3FGL Catalog and our 7-years analysis. Following up with the hypothesis 
of 2FGL J0545.6+6018 being a DM subhalo we interpret its SED in the y-ray band to be purely 
originated by annihilation of WIMPs in the object. We model the differential flux as in Eq. 5.1, 
where ( <J ann v) is the thermally-averaged self-annihilation cross-section, mwiMP is the WIMP mass, 
dNy/dE is the photon yield per annihilation through channel i, IV is the branching ratio for channel 
i, and p represents the density distribution of DM. We assume a value for the thermally-averaged 
self-annihilation cross-section of (o cmn v) = 2.2 x 10 26 cm 3 s 1 [15] and consider the parametriza- 
tions from [16] as coded in the DAMASCO package 4 for the photon yield per annihilation. Thus, 
our model contains two free parameters, mwiMP from the particle physics factor and the astrophys- 
ical factor J. 

0(E,AG) = 0 pp (E)x7(AD) = ^^f^£R'^dEx f p 2 (r(s,Cl))dsdfl . (5.1) 

4 it 2mf M f-[ dt J An, l.o.s. 

5 -v-•' S V ■ / 

Particle physics factor Astrophysical factor 

Four channels are considered: annihi¬ 
lation into bb, W + W~, p + p~, and z + z~, 
with branching ratios B' = 1 in all scenarios. 

While the fit did not converge for the p + p~ 
channel and the goodness of fit for the z + z 
is really poor (the hypothesis of DM annihi¬ 
lating into z particles is rejected at a < 0.01), 
the fits to annihilation into bb, W ] W pro¬ 
vide with decently good fits. The best-fit values 
together with their statistical errors can be found 


Channel 

m WIMP 

logioW 

X 2 /d.o.f. 

p— value 


[GeV] 

[logio(GeV 2 cm- 5 )] 



bb 

78.3 + 11.6 

20.8/q] 

1.02 

0.36 

w+w- 

89.2 ± 13.7 

20.9fjjj 

0.74 

0.48 

1 

D* 

+ 

L* 

18.0 + 0.3 

20.3L8;} 

14.69 

4 x 10- 7 


Table 4: Summary of 2FGL J0545.6+6018 SED fit to 
DM annihilation spectra. We assume (o ann v) = 2.2 x 
10 26 cm 3 s 1 for the estimation of the /-factors. 


for the WIMP mass and the astrophysical factor 
in Table 4. 


4 http://cta.gae.ucm.es/gae/damasco 
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The most favored annihilation channel 
in terms of the goodness-of-fit is IT IT , 
with a Xd 0 f = 0-74, although the fit result is 
in slight tension with the lowest-energy up¬ 
per limit in the SED as showed in Fig. 3. On 
the other hand, annihilation into bb, with a 
Xd 0 f = 14)4, is compatible with the lowest- 
energy upper limit in the SED. We note that 
the later annihilation channel and its associ¬ 
ated WIMP mass of ~ 89 GeV are similar 
to the results presented in [7] for a different 
set of DM subhalo candidates, and in reason¬ 
able agreement with the Fermi -LAT Galactic 
Center excess DM modeling in [17]. The 
estimated 7-factors are one order of mag¬ 
nitude above the nominal 7-factor for the 
most promising ultra-faint dwarf spheroidal 
galaxy in [18], namely Ursa Mayor IP.The 
most constraining limits in this energy range come from the analysis of Fermi -LAT data on dwarf 
spheroidal galaxies [19]. For annihilation into bb and mwrnp = 78 GeV, they found (<7 ann v) < 
1.0 x 10- 25 cm 3 s -1 at 95% cl. Similarly, for annihilation into IT IT and mwrnp = 89 GeV, they 
found (Oannv) < 1.6 x 10 23 cm 3 s 1 at 95% cl. Therefore, although in tension with the latter 
limits, our results are not ruled out by them. 

6. Summary and outlook 

We have presented a search for DM subhalo candidates in the 2FGL Catalog and the results 
from the observations with VERITAS of two promising candidates, 2FGL J0545.6+6018 and 2FGL 
J1115.0—0701, along with a dedicated 7-years LAT data for both sources. 

2FGL J1115.0—0701 flux variability and high Galactic latitude opens the possibility of this 
source being a blazar. We encourage the monitoring of 2FGL J1115.0—0701 HE y-ray activity 
searching for potential flaring events which could drive its VHE y-ray component, depending on 
the source’s redshift, into the sensitivity reach of the current generation of IACTs. 

We interpreted 2FGL J0545.6+6018 SED in terms of DM annihilating via different channels. 
Annihilation into bb and IT IT fit the source’s y-ray emission reasonably well, providing WIMP 
masses of 78.3 GeV, and 89.2 GeV correspondingly. While 2FGL J0545.6+6018 y-ray emission 
has been interpreted in terms of a DM signal, less exotic explanations compatible with its spectral 
shape can be offered, like attributing its emission to a radio-faint y-ray pulsar. Multiwavelength 
follow-up observations are encouraged in order to elucidate the real nature of this enigmatic source. 

It is worth mentioning the impact that the next-generation of IACTs, represented by the 
Cherenkov Telescope Array (CTA, [20]), may have on this kind of searches. CTA is expected 

5 logio(7uMall) = 19.9 ^q 5 for an integration angle of 0.5° (the 7-factor units are GeV 2 cm -5 ). 
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Figure 3: 2FGL J0545.6+6018 spectral energy dis¬ 
tribution. The two best-fit annihilation spectra are 
shown. 
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to improve current-generation IACTs’ sensitivity by a factor of ten which, in conjunction with its 
survey capabilities, will increase the likelihood to serendipitously discover new VHE y-ray sources 
and, potentially, DM subhalos. 
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